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ABSTRACT 

We analysed thermonuclear (type-I) X-ray bursts observed from the low-mass X-ray binary 
4U 1728-34 by RXTE, Chandra and INTEGRAL. We compared the variation in burst energy and 
recurrence times as a function of accretion rate with the predictions of a numerical ignition model 
including a treatment of the heating and cooling in the crust. We found that the measured burst 
ignition column depths are significantly below the theoretically predicted values, regardless of the 
assumed thermal structure of the neutron star interior. While it is possible that the accretion rate 
measured by Chandra is underestimated, due to additional persistent spectral components outside the 
sensitivity band, the required correction factor is typically 3.6 and as high as 6, which is implausible. 
Furthermore, such underestimation is even more unlikely for RXTE and INTEGRAL, which have 
much broader bandpasses. Possible explanations for the observed discrepancy include shear-triggered 
mixing of the accreted helium to larger column depths, resulting in earlier ignition, or the fractional 
covering of the accreted fuel on the neutron star surface. 

Subject headings: X-rays: bursts — X-rays: individual (4U 1728—34) — stars: neutron — X-rays: 
stars 



1. INTRODUCTION 

Thermonuclear (type I) bursts are triggered by unsta- 
ble nuclear burning of the material accreted onto the neu- 
tron star (NS) surf ace in low-ma s s X-ray binary (LM XB) 
syste ms (e.g., iLewin et al.lll995t iStrohmaver fc Bildstenl 
2006). The basic theory of typ e I bursts was outlined 
short ly aft e r their detection (e.g.-IWoosley fc Taamlll976l: 
Joss! 119771: IMaraschi fc Cavalierd Il977t ILamb fc Lambl 



1978). According to these models, the accreted mate- 



rial, which usually consists of hydrogen and helium, ac- 
cumulates as a thin layer on the NS surface (typically on 
time-scales of hours to days) , and when the pressure and 
temperature at its base reach critical values, the fuel will 

ignite and burn unstably un til exhausted. 

Further observations (see Gallowa y etal.l I2008L for a 
revie w) and su bsequent model ing dFuiimoto et al. 
19811 fl987t iFushiki fc Lambl 119871: IBildsten 
199a INaravan fc Hevll 120031: IWooslev et al.l 12001 



Cooper fc Nara van 200(3) showed in detail how the burst 
properties depend on the accretion rate, composition 
of the accreted material (the H/He fraction and CNO 
metallicity), and internal properties of the neutron 
star. For example, for systems accreting mixed H/He, 
the heat generated from hydrogen burning is usually 
a dominant factor for ignition. On the other hand, in 
evolved systems in which little or no H is present but 
the fuel consist mainly of He, the heat required for the 
burst ignition must come entirely from the electron 
captures and pycnonuclear reactions in the NS crust. 
Comparisons of He-bursts with ignition models thus 
offer a powerful probe of the physical conditions in 
the neutron star crust, below t he fuel layer, and the 
cooling processes in the core (jFujimoto et al.l Il987t 
ICumming et alJl2006f >. 

The best-known He-accretor is a low-mass binary sys- 
tem 4U 1820—30, in which the neu tron star orbits it s 
companion once every 11.4 minutes (Ste lla et al.lll987t) . 



Such a tiny orbit cannot accommodate a H-rich com- 
panion, and the mass donor is li kely to be a He-rich 
white dwarf (|King fc Wat son! 1 19861) . 4U 1820-30 is in a 
bursting mode for aroun d 40 days after switching to the 
low state (jChou fc Grind lav 200 j) while during the rest 
of its p» 176-day accretion cycle (|Priedhorskv fc Terrelll 
1984) the source does not exhibit bursts (Stella et al. 
1984 and references therein). In a 20- hour EXOSAT 



observation during the low state, lHaberl et al.l (|1987t ) 
observed nearly regular bursting from 4U 1820—30, de- 
tecting seven bursts with a mean recurrence time of 
3.21 ± 0.04 hours and persistent luminosity of Lx = 
2.8 x 10 37 ergs s _1 between the bursts. 

The b urst ignition cond itions for this source were mod- 
elled by ICumm ing (2003), who compared the predicted 
burst pro perties wit h the m easurements of lHaberl et al.1 
(|1987| ). ICummingl (J2003I ) presented models for pure 
He fuel, but also estimated the effect of adding a 
small amount of hydrogen - 5—35% by mass - as 
predicted by some stella r evolutionary models (e.g., 
iPodsiadlowski et all 120021 ). The amount of energy re- 
leased in pycnonuclear and electron capture reactions 
that escapes from the surface (Qcrust) is a free param- 
eter in this model, which also assumes the time-averaged 
rather than instantaneous accretion rate to set the crust 
temperature profile, because the thermal time in the 
crust is m u ch lon ger than the «176-day accretion cycle. 
ICummingl (|2003f ) found a good agreement between the 
model, which assumes a mixed fuel (10% of hydrogen), 
and the data, providing that Q C rust = 0.1 MeV/nucleon 
(j Brown! 120001 ) and the time-averaged accretion rate is 
rs2 times larger than the measured rate. However, for 
a fuel consisting entirely of He, the required Qcrust was 
0.4 MeV/nucleon and the required time-averaged accre- 
tion rate was 4-5 times larger. Self-consistent results 
were obtained with improved burst ignition models, in 
which the flux flowing outwards was calculated directly 
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from the neutron star crust and core n eutrino emissiv- 
ity a nd the core therm a l cond uctivity (|Cumming et al.l 
I2006D . iCumming et al.l (|2006| ) concluded that, in order 
to produce the bursts separated by «3 hours, the core 
neutrino emissivity must be very inefficient (e.g., supp- 
resed modified Urea process) and the accretion rate must 
be «2 times larger than that inferred from the X-ray lu- 
minosity. 

Although all previous studies of He-bursts have focused 
on 4U 1820—30, the intermittent occurrence of the bursts 
makes triggering of burst observations extremely diffi- 
cult, and only a few recurrence times and corresponding 
accretion rate estimates are available. A much more suit- 
able candidate for such studies is the source 4U 1728—34, 
which consistently exhibits frequent bursts characteris- 
tic of pure He fuel. There is a total of 106 bursts from 
this s ource in the RXTE burst catalogue (|Galloway et al.1 
,2008). The a values («200), short rise times and decay 
time scales suggest a He-rich fuel. The persistent flux 
during RXTE observations was 1 — 7 x 10 -9 erg cm~ 2 
s _1 while the burst recurrence times were on average ~4 
hours. Evidence for a short orbital period of 10.77 min 
has been detected recently in the analysis of Chandra 
observations (Galloway et al. 2010; in prep.), supporting 
the long-suspected identification of 4U 1728—34 as an 
ultracompact LMXB. The source is probably accreting 
pure He from its evolved companion and is practically a 
twin of 4U 1820—30, except for the much more frequent 
and reliable bursting. 4U 1728—34 has also been ob- 



served extensively by I N TEGRAL (IFalanga et al.l 120061: 
iChelovekov et al.)l2006D . IFalanga et all (|2006D detected 



36 type I bursts during the transition from hard to soft 
state, where the source luminosity increased from 2—12% 
of the Eddington luminosity. 

In this paper we present analysis of new Chandra ob- 
servations of 4U 1728—34. We measured the recurrence 
times and corresponding accretion rates of the 25 bursts 
detected during a 240-ks HETGS exposure. The detailed 
spectral analysis and the detection of the orbital period 
and radius expansion bursts are reported in the compan- 
ion paper (Galloway et al. 2010), while we present the 
comparison of the observed burst properties with a new 
ignition model. In Section 2 we present the data and 
describe our analysis of the Chandra observations. The 
measured burst properties, which include the accretion 
rates (estimated from the persistent flux), burst fluences, 
a values and burst recurrence times are presented in Sec- 
tion 3. The new burst ignition model is described in Sec- 
tion 4, while Section 5 shows the comparison of the model 
and data. Possible explanations of the discrepancy be- 
tween the data and our ignition model are discussed in 
Section 6. Finally, our conclusions are summarized in 
Section 7. 

2. OBSERVATIONS AND ANALYSIS 

2.1. Chandra observations 

4U 1728-34 was observed between 2006 July 17-23 
with the HETGS abroad Chandra. The observations 
6568, 6567 and 7371 were made in the continuous clock- 
ing (CC) mode with exposures of 49.5, 151.8 and 39.7 
ks, respectively. A total of 25 type I bursts, separated by 
1.8—3.9 hours, were detected in these observations (for 
the light curve and more details on the Chandra obser- 



vations see Galloway et al. 2010). Galloway et al. (2010) 
also report on the detection of a period of 10.77 min in 
the low-energy persistent intensity (which they interpret 
as arising from orbital modulation) and describe a search 
for line emission and photoionization edges in the persis- 
tent spectra between the four radius expansion bursts 
detected in observation 6568. In this paper we focus on 
measuring the pre-burst accretion rates and correspond- 
ing burst recurrence times, and their comparison with 
our new ignition models. 

The first order HEG and MEG persistent spectra were 
extracted using time intervals starting 150 seconds af- 
ter the previous burst, and ending 50 seconds before 
the burst peak, to exclude any possible burst emis- 
sion. These spectra, which contained approximately 
2000 to 3000 total counts, were binned to a minimum 
of 100 counts per bin, and fitted simultaneously for 
each interval. The combination of the CC mode (se- 
lected to minimize photon pileup) and grating observa- 
tions is non-standard, making the correct background 
subtraction possible only for the brightest sources (e.g., 
iMarshall et al.l 120011 ) . Hence, we restricted our spectral 
fitting to the 1.5—6 keV band, in which the contribution 
from the background is minimal. The response matri- 
ces were produced by the CIAO task mkgrmf, while the 
task fullgarfs was used to produce the auxiliary files, 
which were used for fitting both the persistent and time- 
resolved burst spectra. 

Since the Chandra energy band is relatively narrow, to 
help us select the appropriate spectral model, we exam- 
ined spectral fitting results of previous observations made 
in significantly broader bands. The source was in a soft 
state during t he broadband (0.1 — 1 00 keV) BeppoSAX 
observation by Di S alvo et aD (|2000T i . The best- fit model 
consisted of two components, a 2-keV blackbody (with an 
emitting region comparable to the expected radius of the 
neutron star), and a Comptonized component (comptt 
model in xspec) with seed photon temperature of w 1.5 
keV, ele ctron temperature of ~ 10 keV and optical depth 
of «5. iDi Salvo et al.1 (|2000f ) also detected two broad 
emission lines at ~6.7 ke V and sal .6 keV, pr obably emit- 
ted in the ionized corona. iD'Af et aLl (J2006I ) reported on 
the analysis of the simultaneous RXTE and Chandra ob- 
servations of 4U 1728-34. The broadband (1.2-35 keV) 
spectrum was best described by a blackbody {kT «0.6 
keV) and a Comptonized component (fcTo wl.5 keV; 
kT e «7 keV; r ~5), but no emission lines were detected. 
Instead, ID'Af et all (|2006l ) have found absorption edges 
at «7 keV and «9 keV associated with Fei and Fexxv. 
During several INTEGRAL observations (in the 3-200 
keV band), the source was found undergoing the tran- 
sition from the intermediate/hard to soft state with the 
electron temperatures of the Comptonizing plasma de- 
creasing from «35 keV to 3 keV (jFalanga et al. 2006). 

Following these results, we selected a single- 
component 1 absorbed comptt model in XSPEC to fit 
the HEG and MEG data. We then used the spectral 
results from previous observations to fix some of the 

1 A blackbody component may also be present when the source 
is in the soft state, in which case the estimated bolometric flux 
would be reduced (by at most fs 20% for our data). By including 
a blackbody component the contribution from the comptt com- 
ponent is reduced, and the latter contributes more to the total 
bolometric flux, particularly at high energies. 
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model parameters as the rather limited Chandra band 
was not sufficient to constrain them. First, we adjusted 
and fixed the seed photon energy to 0.4 keV for the first 
four bursts detected in the observation 6568, 0.5 keV 
for the 18 bursts in 6567 and 0.3 keV in the last three 
bursts detected in the observation 7371. These values 
were selected so that the corresponding absorption col- 
umn was consistent with the best-fit value found from the 
burst spectra (N H — 2.29 x 10 22 cm" 2 , see below), al- 
though we actually found that the measured fluxes were 
not very sensitive to this parameter. Since the Chandra 
band was well below the spectral turnover (wlO-20 keV, 
for kT e ?«3 keV), we could select and fix a wide range 
of electron temperatures and get statistically acceptable 
fits. We selected the temperature of 35 keV to represent 
a hard spectral state, and then repeated the fitting with 
a low (3 keV) plasma temperature typical of a soft state. 
Since the narrow Chandra band was not sufficient to dis- 
criminate between these spectral states, we used both 
models to estimate the hard and soft-state persistent X- 
ray fluxes for each burst, from which we calculated the 
corresponding accretion rates. 

To measure the neutral hydrogen column density, we 
first produced the peak spectra of all bursts using only 
the data around the peak (between 0.5 s before the 
peak and 1.5 s after the peak) and combined these spec- 
tra for each observation. Since Galloway et al. (2010) 
showed that the burst profiles in each observation have 
approximately the same shape (see their Fig. 3), we 
combined the peak spectra of the four bursts from the 
observation 6568, eighteen bursts from 6567 and three 
bursts from 7273. We subtracted the background (per- 
sistent) emission from these combined spectra and fit- 
ted them with absorbed blackbody models. We ob- 
tained good fits for each of the three observations with 



an N H = 2.30^ 



x 10 cm in the first observa- 
tion, N H = 2.27t°;44 x 1Q22 cm ~ 2 in tlie second, and 
N H = 2.3ll 1 i3 x 10 22 cm" 2 during the last three bursts. 
We then calculated the weighted mean and variance of 
iV H = 2.29±°; 2 ,f x 10 22 cm" 2 . 

Finally, we produced time-resolved spectra for each in- 
dividual burst by applying an adaptive time-binning, i.e., 
by accumulating the counts until a minimum of approx- 
imately 800 counts (total HEG and MEG first order) in 
each bin. We then fitted the HEG and MEG spectra 
of each bin with an absorbed blackbody model (with the 
Nh fixed at the mean value of 2.29 x 10 22 cm -2 ) and used 
the model parameters to calculate the corresponding 
bolometric flux. The fluence of each burst was then cal- 
culated by summing the measured fluxes over the burst. 
As one example, in Fig. [T] we show the best-fit blackbody 
model parameters and fluxes for time-resolved spectra of 
two bursts from our sample. Since the absorption col- 
umn measurements from persistent spectra in previous 
observations range from 1.6 — 2.8 x 10 22 cm" 2 , we esti- 
mate the maximum possible systematic error of 20 — 25%. 
We also note that for the RXT E burs t s sele cted for our 
analysis (Section I2.2p Gall oway et al.l ((2008) performed 
the spectral fitting with the absorption column fixed at 
the mean value for each burst. For the 11 RXTE bursts 
included here, this value is in the range 0.6 — 2.4 x 10 22 
cm" 2 . 
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Fig. 1. — The fitted blackbody temperature, radius (assuming 
a distance of 5.2 kpc) and bolometric flux measured during the 
bursts #1 (top) and #5 (bottom) in the Chandra data. The error 
bars indicate the 1 a uncertainties. Each time-bin contains ss800 
cts. The burst peak occurs 1—1.5 sec from the start of the burst 
(t=0 sec). 



2.2. RXTE and INTEGRAL bursts 

From the total of 106 burst s observed from t his so urce 
by RXTE and re-analysed bv iGallowav et al.l (|2008[ ). we 
selected 11 bursts for which the recurrence times could 
be reliably measured. We used the catalogued values of 
the persistent X-ray luminosities (measured in the 2.5- 
25 keV band and multiplied by the bolometric correction 
factor of 1.38; Galloway et al. 2008), burst fluences and 
a parameters. 

4U 1728-34 was in the field of view of JEM-X and 
ISGRI cameras during several thousand short INTE- 
GRAL observations of the Galactic centre region made 
in 2003 and 2004. A total of 36 bursts detected 
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Fig. 2. — Top: The burst fluence versus the accretion rate mea- 
sured from the persistent X-ray luminosity and expressed in units 
of the Eddington luminosity, assuming the accreted material to be 
pure helium. Bottom: The burst fluence versus the recurrence 
time. 

in the co mbined INTEGRAL observations were cata- 
logued by Fal anga et al.l (|2006f> . while a complete anal- 
ysis of all availab l e INT EGRAL IBIS observations by 
IChelovekov et all (|2006[ ) revealed 61 bursts from the 
source. The burst recurrence times, however, could only 
be accurately measured if the source was in the JEM- 
X field of view continuously between two subsequent 
bursts (M. Falanga, private communication), which we 
could verify for 5 bursts. These 5 bursts, which we 
include in our an alysis, are listed in t h e INTEGRAL 
burst catalogue by iFalanga et al.l (|2006l ). Falang a et al.l 
(|2006f ) measured the X-ray luminosities (in the 3—200 
keV band) and used them to estimate the accretion rates. 
The burst fLuences were measured from the time-resolved 
burst spectra with a typical bin size of 1—2 seconds. 



3. MEASURED BURST PROPERTIES 
3.1. The persistent luminosity and accretion rates 

Using the best-fit models of the persistent pre-burst 
spectra, we measured the absorbed persistent fluxes in 
the 1.5—6 keV band to be in the range 1.16 — 2.30 x 10~ 9 
ergs cm -2 s _1 , and estimated the corresponding intrinsic 
bolometric luminosities 2 (1.3 — 2.3 x 10 37 ergs s _1 for the 
fcT c =3 keV, and 4.1 - 6.4 x 10 37 ergs s" 1 for the fcT c =35 
keV), from which we calculated the pre-burst mass ac- 
cretion rates. The range of the estimated bolometric lu- 
minosities is approximately 10% smaller than that of the 
fluxes in the 1.5-6 keV band, due to the different choice 
of kTo for the three observations (0.3, 0.4 and 0.5 keV for 
observation 7371, 6568 and 6567, respectively; see Sec- 
tion 2.1). Simple phenomenological (power-law) fits indi- 
cate that the persistent spectra actually become slightly 
harder with increasing flux, suggesting that the range 
in bolometric flux should be larger rather than smaller. 
However, these effects are relatively minor compared to 
the uncertainty in the electron temperature T e . 

The accretion rates calculated from the persistent 
flux, assuming Comptonization in the high temperature 
plasma (kT e «35 keV), would imply almost one order of 
magnitude larger a parameters (Section 13. 2[) than usu- 
ally measured from this source. In addition, since the 
relatively high persistent flux during the Chandra obser- 
vations suggests a soft spectral state, we conclude that 
the accretion rates estimated from the 3-keV X-ray lu- 
minosities are more realistic, and we will assume these 
values throughout. 

The accretion rates are expressed in units of the Ed- 
dington limit. The Eddington accretion rate for pure he- 
lium fuel is assumed to be 1.78 x 10 18 g s _1 in the distant 
observer's frame, which corresponds t o the Eddington 
lumin osity of 3.5 x 10 38 ergs s^ 1 (e.g., iNaravan fc Hevll 
l2003f h Although the ultra-compact systems are thought 
to accrete almost pure helium (i.e., the fraction of the 
accreted hydrogen X — 0), some stellar evolution mod- 
els p redict that the hydrogen f raction can be as high as 
20% (jPodsiadlowski et"aT1l2002h . To account for this pos- 
sibility, wc compared all our calculations with the case 
Xq = 0.2, where we used the corresponding Eddington 
accretion rate of 1.48 x 10 18 g s^ 1 (L Ed d = 2.9 x 10 38 
ergs s" 1 ; INaravan fc Hevll 120031) . 

The persistent accretion rate during the first two Chan- 
dra observations was nearly constant at around 0.066 of 
the Eddington limit, decreasing to «0.038 during the last 
pointing, when the last three bursts occurred. The per- 
sistent accretion rate during the five bursts detected with 
INTEGRAL and eleven RXTE bursts varied between ap- 
proximately 0.02 and 0.08 of the Eddington accretion 
rate (M E dd.)- 

3.2. Burst fluences and a parameters 

Fig. [2] shows the burst fluences, measured from the 
time-resolved burst spectra (Section [2|), and plotted 
against the burst accretion rate (left) and recurrence time 
(right). The bursts are weak (1 — 8 x 10 39 erg) and we 
see a sharp increase of the burst energy with burst re- 
currence time, but only for short recurrence times (~ 2 

2 assuming a distance of 5.2 kpc measure d from the photospher ic 
radius expansion bursts in RXTE data by Galloway et al. (2008). 
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hours) . 

The parameter a is usually used to compare the gravi- 
tational energy (Q gra v) produced by accretion and the 
energy released by nuclear burning (Q nuc ) during the 
burst (a = zc 2 /Q nuc ] where z is the surfa ce gravitational 
redsh ift and c is the speed of light; e.g., iFuiimoto et al.l 
1987). The nuclear energy depends on the fuel compo- 
sition and is calculated as Q nuc = 1.6 + 4(X) MeV per 
nucleon, where the hydrogen fraction in the fuel, (X), 
is averaged over the burning layer. The formula gives 
an energy of 4.4 MeV per nucleon for the Solar abun- 
dance of hydrogen ((X) — 0.7), which implies 3 a p» 60, 
while higher values suggest a He-rich fuel. We have 
calculated the a parameters for the Chandra bursts as 
a = F p x At/Eb, where F p is the persistent flux inte- 
grated over the burst recurrence time At, and E\, is the 
burst fluence, and found it to be in the range 105—342 
(the mean value 189), which implies a He-rich fuel. Con- 
sistent a values are catalogued for the RXTE (91—314) 
and INTEGRAL bursts (160-240) with the mean value 
of 180 for the whole sample. 

In addition to the a parameter, from the measured 
burst fluences (E\>) we calculated the total mass burnt 
during the burst (Mb = Eb/Q nuc ; where Mb is mea- 
sured in the distant observer's frame). This allows us to 
compare the consumed mass during the burst directly to 
the mass accreted since the preceding burst. According 
to the thin-shell instability models, the accreted mass 
is burnt completely during the bursts, unless there are 
some energy "leaks", for example, due to stable hydro- 
gen burni ng before the burst o r other sources of energy 
loss (e.g..fFuiimot o et aLlll987l and references therein). 

We calculated Mb for the pure helium fuel, and re- 
peated the calculation for Xq = 0.2. For the Solar 
metallicity, this fraction of hydrogen will be completely 
burnt in ss3 hours, via the hot CNO cycle (e.g.. lCummind 
2004) . Since most of the measured burst recurrence times 
are between 3 and 5 hours (Fig. [2]), no hydrogen is left 
to burn unstably and the Q n uc does not increase signifi- 
cantly for these bursts. However, if the metallicity is, for 
example, ten times lower (Z = 0.0012), the time needed 
to burn the accreted hydrogen is almost 32 hours, which 
means that most of the accreted hydrogen will be left to 
burn unstably during the burst. As a consequence, Q n uc 
will increase up to ~ 2.3 MeV per nucleon for bursts with 
short recurrence times. 

Fig. [3] shows the ratios of the accreted and burnt mass 
and corresponding a parameters versus the measured ac- 
cretion rate, assuming the fuel consists of pure helium. 
The accreted mass is on average 1.4 times larger than the 
total mass of fuel burnt during the burst. However, there 
is a large scatter around this value, with the fraction of 
the exhausted fuel of «50% or less during several bursts. 
We also notice several bursts with the total burnt mass 
that appears to exceed the accreted mass (i.e., the mass 
ratio < 1). Mass ratios greater than «2 are measured for 
three bursts, which also have large a parameters («300). 
A small a value of «90 was measured for the two bursts 
with the corresponding mass ratio of «0.7. 



The 



approximate 



formula, 



-i 



44 MeV/nuc (M/lAM sun )(RJ 10 km ) " 1 (Q nuc /4.4 MeV/nucl) " 

is given by Galloway ct al. (2008), who omitted the redshift 
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Fig. 3. — The ratio of the accreted mass and the mass consumed 
during the burst plotted as a function of the persistent accretion 
rate (in units of the Eddington limit, assuming the accreted mate- 
rial to be pure helium) for Chandra (diamonds), RXTE (asterisks) 
and INTEGRAL data (squares). The corresponding a parameter is 
plotted on the right Y-axis. The line shows the a values predicted 
by our model (see Section 5). 

4. IGNITION MODELS 

We use the general-relativistic gl obal linear stability 
analysis of lCooper fc Narayanl (J2005) to theoretically de- 
termine the type I X-ray burst recurrence time and a 
value as a function of accretion rate. We assume steady 
spherical accretion onto a neutron star of gravitational 
mass M = 1.4 M© and areal radius R = 10.4 km at 
a rate M, where M is the rest mass accreted per unit 
time as measured by an observer at infinity. We describe 
the accreted matter's composition by the hydrogen mass 
fraction X, metal mass fraction Z , and helium mass frac- 
tion Y = 1 - X - Z. 

We make the followi n g two modific ations to the model 
of iCqoper fc Naravanl (|2005l ). (1) iCooper fe Naravanl 
(2005) assumed the energy generated by electron cap- 
tures, neutron emissions, and pycnonuclear reactions in 
the crust was dis t ribute d uniformly. We now follow 
IHaensel fe Zdunikl (J2008D and distribute the energy ac- 
cording to their Table A. 3. (2) We updated the crust 
therm al co nductivity according to IShternin fc Yakovlevl 
(|200H and lChueunov fc Haensell (|2007l) . 

The accretion rate, accreted matter composition, and 
interior temperature profile together determine the recur- 
rence time. The temperature profile is a strong function 
of the crust's thermal conductivity and core's neutrino 
emissivity, both of which are poorly constrained. The 
thermal conductivity is a decreasing f unction of the im- 
purity para meter Q 



mm 



(Z) 2 (lltoh fc Kohvama 
1993: Daligaul t fc Guptall2009l ). Although iSchatz et al. 
(1999J) found Q 



imp 



correction factor 1 - 



1.3. 



100 from steady-state nucle- 
osynthesis calculations, subsequent investigations sug- 
gest Qimrj should be much smaller (ISchatz et ail 120031: 
i Wooslev et"aT] 120041: iKoike et al.lj2004t iHorowitz et al. 
2007, 2009: Shternin et al.1 120071: L Brown & Gumming 
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2009). Also, the core neutrino emissivity, and thereby 
the core cooling rate, depends on the unknown ul- 
tradense matter equation of state ( for reviews, see 
lYakovlev fc Pethickl l200l iPage et al.l l2006h . Unless 
noted otherwise, we set Qi mp = 100 throughout the crust 
and adopt a slow, suppressed core cooling model in which 
nucleon-nucleon bremsstrahlung proce sses dominate th e 
neutrino emission (see, e.g., Table 1 of IPage et aLl [2006): 
these parameters are realistic limits on the crust ther- 
mal conductivity and core neutrino emissivity and are 
set to maximize the burst ignition region temperature 
and thereby minimize the burst recurrence time, for the 
closest correspondence with observations. 

5. COMPARISON WITH OBSERVATIONS 
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Fig. 4. — Left: The burst recurrence time (top), the correspond- 
ing accreted column depth (calculated by multiplying the local 
accretion rate by the burst recurrence time; middle) and burst ig- 
nition column depth (calculated from the burst fluences; bottom) 
plotted as a function of the accretion rate. The accretion rate is 
in units of the Eddington accretion rate assuming pure helium fuel 
(Xq = 0). The line shows the model: a slow, suppressed neutrino 



emission process in the core and Qi n 



100, assuming the ac- 



creted material to be pure helium. The error bars represent the 
data: Chandra (diamonds), RXTE (asterisks) and INTEGRAL 
(squares). The mass accretion rates are inferred from the persis- 
tent fluxes measured in the energy bands of the three instruments 
(see text for details). Right: The same plot for the accreted ma- 
terial containing 20% of hydrogen. The corresponding Eddington 
accretion rate is reduced. The two models shown are again for the 



slow, suppressed cooling neutrino process with Q; n 



100, but 



assuming the accreted material has Xq = 0.2 and a Solar metallic- 
ity (full line), and a metallicity ten times less than the Solar value 
(dashed line). 
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Fig. @] (left) shows the model and the observed burst 
recurrence time and ignition column depths (calculated 
from the accreted mass and also from the mass consumed 
during the burst) plotted versus the mass accretion rate 
for Chandra, RXTE and INTEGRAL data assuming the 
accreted material to be pure helium. Although the se- 
lected model maximizes the outward heat flux that ig- 
nites the bursts, the comparison of the observed and the- 
oretical burst recurrence times shows a large discrepancy. 
For example, for a typical accretion rate of ~ 0.05MEdd 
we observe bursts separated by 3—4 hours, while the 
model predicts almost 15 hours. To match the observed 
burst recurrence times each of our estimated accretion 
rates would have to be multiplied by a correction factor 
in the range 2.5-6 (mean 3.6). Such a large discrep- 
ancy may suggest that the region at the base of the igni- 
tion column is significantly hotter than predicted by the 
model, which means either that the cooling processes in 
the crust and core are slower than the current theoretical 
lower limits, or that there is an additional heat source. 

The additional heat may be supplied by burning of hy- 
drogen, if the accreted material is a mixture of hydrogen 
and helium. Since the stellar evolutionary models pre- 
dict that a small fraction of hydrogen might be present 
in evolved, ultracompact binary systems, we produced a 
model assuming the same core and crust cooling proper- 
ties, but with the accreted material that includes around 
20% of hydrogen (Xq = 0.2), which is a reasonable up- 
per limit for an ultracompact binary. Since the hydrogen 
burning is affected by the CNO fraction, we produced 
two models, one with the Solar CNO fraction of 0.012, 
and the other with the CNO fraction ten times lower 
(Z=0.0012). We also re-calculated the Eddington accre- 
tion rate limit for the new fuel composition, and the new 
models are shown in FigQ] (right). Although the discrep- 
ancy between the model and data was slightly reduced, 
adding 20% of hydrogen was not sufficient to match the 
theoretical and the observed burst recurrence times. The 
model with the low CNO fraction predicts shorter burst 
recurrence times. This is expected, since the assumed 
hydrogen fraction of Xq = 0.2 needs more than 30 hours 
to burn stably when the metallicity is 10 times lower 
than the Solar CNO fraction (see Section |3T2]) . Hence, 
the hydrogen will burn stably for longer time and ignite 
later, at a larger depth and hence closer to the helium 
burning region 4 . 

To estimate the "missing heat" we produced several 
models, which further increase the outward heat flux. 
These models also assume the slow cooling and low con- 
ductivity of the neutron star core and crust as before, 
but include additional heating. As our first attempt, we 
increased the total heat produced deep in the crust by 
electron capture and pycnonucear rea ctions from the cur- 
rent theoretical estimate of 1.9 MeV (Hacnse j fe Zdunikl 
12008ft to 3.0 MeV per nucleon. Although this slightly 
reduced the burst recurrence times, it was not enough to 
match the observed values. 

Recently, iBrown fc Cummingl ()2009h investigated the 
cooling rates of several quasi-persistent transients dur- 
ing the quiescence. They found that, in order to fit 

4 This is not the case, however for very large mass accretion 
rates, because the ignition of helium occurs at much lower column 
depths 



the cooling light curves of these sources following an 
outburst, an inward-directed heat flux from the shal- 
low outer cr ust was required. For sev eral sources in 
their sample, IBrown fc Cummingl (|2009ft estimated the 
strength of this shallow heat source to be in the range of 
0.7—1.1 MeV per nucleon. Althoug h the nature of this 
shallo w heat source is not yet clear, IBrown fc Cummingl 
(2009) suggest that its strength probably depends on the 
accretion rate. 

Motivated by this result, we included an additional 
shallow heat source in our model. For a strength of 0.7 
MeV per nucleon, the additional heat source's effect was 
larger than that from the stronger heat produced deeper 
in the crust: it reduced the burst recurrence time from 
w 10 hours to w 6.5 hours (for M « 0.07AfEdd), approxi- 
mately two times longer than observed. However, regard- 
less of the heat source's strength, the predicted ignition 
depth was always > 1 x 10 8 g cm -2 , which is well above 
the observed values (Fig. [3J. Indeed, we find that helium 
burns stably for a sufficiently large assumed strength of 
the shallow heating. The reason is the following (see, 
e.g.. lPaczynsk3ll983t fFushiki fc Lamblll987l ). Increasing 
the heat source's strength raises the burning layer's tem- 
perature, which decreases the column depth a t which he - 
lium burning is thermally unstable (e.g., Bildstcn 1998). 
However, helium must survive to that ignition depth; 
increasing the temperature raises the burning rate and 
hence decreases the helium depletion depth. If the former 
exceeds the latter, helium burns stably. 

6. DISCUSSION 

6.1. Comparison with the ignition model 

The most striking feature of the thermonuclear bursts 
exhibited by 4U 1728—34, apart from their consistently 
short rise and decay times, is the short recurrence times 
compared to the predictions of the ignition model. The 
simplest way to explain this discrepancy is that the mea- 
sured accretion rates are underestimated, for example if 
additional soft or hard spectral components are present 
outside the instrumental bands 5 . This is certainly plau- 
sible for the very narrow Chandra band, but additional 
hard components are not possible for the RXTE or IN- 
TEGRAL measurements, unless they occur at very high 
energy (> 100 keV). Such features have not been de- 
tected in LMXBs. Low energy spectral components ap- 
pear more likely, and could easily be overlooked thanks 
to the combination of low-energy sensitivity limits for 
RXTE and INTEGRAL of 2-3 keV, the restriction of 
the Chandra spectral fits to > 1.5 keV, and the strong 
absorption towards 4U 1728—34. However, we note 
that similar disc r epancies were found f or 4U 1820—30 
(| Cummingl 120031 : iCumming et al.l 120061 ) , which is also 
well-studied and has low line-of-sight absor ption (Nh — 
a few times 10 20 cm" 2 ; iMiller et all [2009ft . Studies of 
the burst recurrence time - persistent flux relationship 
in GS 1826—24 indicate that soft spectral components 

5 Another possibility is that the assumed distance of 5.2 
kpc, measured from the photospheric radius expansion bursts by 
Galloway et al. (2008), is underestimated. Although no other in- 
dependent estimates of the distance are available for this source, 
we note that the assumed Eddin gton limit is clos e to th e empirical 
critical luminosity measured by Kuulkcrs et al. (2003). In addi- 
tion, to match the model, the distance to 4U 1728—34 would have 
to be underestimated by a factor of 2—3. 
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contributing up to 40% of the broadband pe rsistent flux 
can be present (e.g.. iThompson et al.l 120081) . However, 
we note that for 4U 1728—34 this additional flux would 
have to be typically 3.6 and up to 6 times larger than our 
measured bolometric fluxes (see Section 5), which seems 
implausible. Furthermore, unless the burst flux was also 
underestimated by a similar degree, such inflated lumi- 
nosities would lead to intrinsic a values that were higher 
by the same factor, presenting additional difficulties. 

An obvious way to increase the burst ignition columns, 
without increasing the global accretion rate, would be 
to restrict the neutron star area on which the nu- 
clear burning takes place. The idea that strong mag- 
netic fields could confine the accreted matter to a frac- 
tion o f t he NS surface was fir st suggested by Uoss fc Lil 
(|1980f) . lAvasli fc Jossl (11982ft later expanded on this 
idea (see also Bilds ten fc Bro wn 1997; Br own fc Bildstea 



1998) although the observationally-inferred dipole mag- 
netic fields of w 10 8 - 10 10 G (e.g.. ILewiti et~a!l Il99l 
iPsaltis fc Lamblll999f ) seem to be much lower than those 
required for the magnet i c conf inement. 

Recently, lLamb et~aH (|2009t ) proposed a model to ex- 
plain the observed properties of the accretion-powered 
oscillations detected in the accreting millisecond X-ray 
pulsars (AMXPs). According to this model, the accreted 
material is channeled to the regions (spots) on the NS 
surface near the magnetic poles, which are close to the 
rotation axis (i. e., the magnetic an d rotational axes are 
nearly aligned). lLamb e t al. (2009) argue that the mech- 
anism that produces the nearly aligned axes could signif- 
icantly reduce the dipole component of the stellar mag- 
netic field, without reducing the the field's total strength, 
which could be as high as 10 11 — 10 12 G (i.e., strong 
enough to confine the accreted material). Furthermore, 
the changes in M and the inner accretion disk structure 
alter the size and location of the X-ray emitting spots, 
which explains the low (sometimes undetectable) pulse 
amplitudes, nearly sinusoidal waveforms and other prop- 
erties of the accretion-powered oscillations from AMXPs. 

The effect of the accreted mass confinement on the 
burst ignition needs to be modelled in detail taking into 
account various effects (e.g., the non-radial thermal diffu- 
sion is zero in the spherically-symmetri c case but nonzero 
in th e confinement case; see, e.g., iPavne fc Melatosl 
2006). Here, we make a simple assumption that the igni- 
tion condition will not change significantly if the accre- 
tion is restricted to only part of the neutron star surface, 
and investigate if the simple area scaling would be plau- 
sible for 4U 1728-34. 

Fig. [5] (left) shows the radius of the restricted area, 
which would be required so that the recurrence time (i.e., 
the accreted column) would match our model shown on 
the left panels in Fig. [4] The radius of the restricted area 
changes from approximately 4.5—6 km, which seem to 
be consistent with the measurements of the black-body 
emitting regions from the time-rcsolvcd burst spectr a 
(e.g., see Fig. [TJ see also Fig. 5 in lGalloway et al.ll2008h . 
Fig. \5\ also suggests that the size of the restricted area 
may be correlated with the global accretio n rate. This 
is pre dicted by the moving spot model of lLamb et al.l 
(2009), suggesting that, perhaps, the nuclear burning in 
4U 1728-34 is also concentrated to moving regions near 
the spin axis. 



In general, a likely consequence of confinement of the 
nuclear burning is that strong oscillations should be ob- 
served during every burst. However, burst oscillations 
were detected in only six out of 1 1 RXTE bursts include d 
in our sample (for more detail see lGallowav et aLi r2008'). 
The absence of burst oscillations during some bursts may 
also be explained if the burst oscillations also arise from 
burning restricted to a region centered close to the spin 
axis and the latitude of this region varies slightly from 
burst to burst. Thus, the absence of strong burst oscilla- 
tions during some of the bursts observed by RXTE from 
4U 1728—34 does not necessarily rule out fuel confine- 
ment. 

In addition to confining the accreted fuel to a smaller 
area, the area scaling would also work if we assume a 
neutron star with a smaller radius. In addition to ac- 
cumulating the fuel over a smaller area, we also expect 
the additional effects (e.g., the higher gravity makes the 
pressure higher for a given column depth, and the higher 
pressure usually gives a higher burning rate, etc.). How- 
ever, based on our calculation for the neutron star radius 
of 6 km, these additional effects are not significant. 

6.2. The observed variations in the a parameter 

Here we address the problem of the large scatter of 
the a values measured for 4U 1728-34 (see Figs. [3] and 
0|). The right panel of Fig. [5] shows the ignition column 
scatter after the area scaling. The model column and 
the column calculated from the measured accretion rate 
coincide, after scaling by the appropriate area (shown on 
the left panel). However, the ignition column calculated 
from the observed burst fluences and scaled by the same 
area shows a scatter. 

Similar large a de viations are ob served in other 
X-ray bursters (e.g., iFuiimotol 119881 . and references 
therein), and a mode l to explain it was suggested by 
Fuii moto et al.1 (|1987| ). According to their model, the 
observed large a deviations from burst to burst can be 
explained by a fuel buffer in the envelope. If the size of 
the fuel buffer is comparable to the total energy of the 
burst (which could easily be the case for the weak bursts 
observed from 4U 1728-34) the a deviations c o uld b e 
large. However, as suggested bv iFuiimoto et all (|1987h . 
the influence of the buffer can be reduced statistically 
if a large number of bursts is observed. If this mecha- 
nism is responsible for the observed a scatter, for con- 
secutive bursts we would expect a pattern of alternating 
high and low a values following each other. Although 
our sample includes 13 consecutive Chandra bursts, the 
large a uncertainties may have prevented us to observe 
this behavior. The a uncertainties of the RXTE bursts 
are significantly smaller, but our sample does not include 
more than three consecutive RXTE bursts. 

The mean value of the a parameter measured from the 
38 bursts in our sample is 180, and we do not see any 
variations of the mean value with the accretion rate, at 
least in our observed range of the accretion rates from 
0.01 — 0.07 MEdd- This is close to the value predicted 
by our model (180-190; see Fig. [3]), which seems to in- 
clude a small "leak" (less than 10%) due to some stable 
burning before each burst. Such a close agreement be- 
tween the average and theoretical a values implies that 
the accretion disk inclination angle with respect to the 
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line of sight is approximately 60 degrees (or slightly less 
for X = (0.2)), as predic ted in the mode l of angular dis- 
tribution of radiation bv lFuiimotol ((19881 see their Fig. 2 
and formula 6). The model suggests that the emission 
is enhanced in the direction perpendicular to the accre- 
tion disk plane, due to scattering in the inner disk, and 
the degree of anisotropy is larger for the persistent emis- 
sion. The angle of «60 degrees is consistent with the fact 
that the eclipses were never observed from this source. 
iFuiimoto et al.l (J1987I ) also argue that the a deviations 
are largest for the photospheric radius expansion bursts 
(PRE) because part of the outer envelope can be ejected, 
increasing the a significantly. Our sample includes four 
possible PRE bursts detected by Chandra (see Galloway 
et al 2010) and also two P RE bursts in the RXTE sam- 
ple (jGallowav et al.ll2008D . but their a parameters are 
not exceptional. This may suggest that the mass loss due 
to the envelope ejection during PRE bu rsts is negligible 
comp ared to the fuel buffer size (e.g., see [Weinberg et ail 
120061) . 

The same mechanism that pushes the envelope buffer 
closer to the ignition zone (i.e., the turbulent mixing and 
dissipative heating associated with thermonuclear insta- 
bilities, cause d by inflow of angul ar momentum with 
accreted gas; IFuiimoto et al.l [1987) could also explain 
weak, frequent bursts ("premature ign ition"). A simi- 
lar ide a was also proposed recently by iPiro fc Bildstenl 
(2007). The turbulent mixing model alone could, per- 
haps, explain the observed burst properties, without in- 
troducing the accreti on confinement, but alt hough its 
spin is slow (363 Hz; iStrohmaver et al.lll996l ). we sus- 
pect that the observed accretion rates of 4U 1728—34 
might be too small for this mechanism to work. Ac- 
cording to this model, the mixing is most effective for 
large accretion rates a nd small spin rates. Figure 15 of 
IPiro fc Bildstenl (|2007h suggests that for the slow spin 
of 4U 1728—34, the required accretion rate to overcome 
the buoyancy barrier would be > O.lAfEdd- However, the 
critical accretion rate depends sensitively on the condi- 
tions in the fuel l ayer and the strength of mixing (e.g., 
see equation 59 of lPiro fc Bildstenl I2007T ). 

It is worth noting that our sample is highly selective for 
bursts with short recurrence times; the bursts observed 
with Chandra all have unusually short recurrence times, 
and for RXTE and INTEGRAL we only quote measure- 
ments for pairs of bursts for which there are no (or few) 
data gaps. If we relax this requirement, there is evidence 
for much higher (up to rs 1000) values of a in the liter- 
ature for 4U 1728—34 and other systems. iFalanga et al.l 
(2006) quote a values up to almost 800, but these are 
upper limits only, since one or m ore data gaps occurred 
in between the bursts. However, I Galloway et al.1 (2008) 
also found that the average a- value for a group of bursts 
with short recurrence times (including 4U 1728—34) was 
> 1000. Such high a values likely cannot arise from in- 
complete burning leaving behind a residual fuel buffer, 
since the required size of the buffer would be many 
times the amount of fuel burnt in the burst. Instead, 
such inefficient burning must arise from a more substan- 
tial energy leak, perhaps steady He-burning, as has also 
been suggested to play a role in the decrease in burst 
rate for most bursters at an accretion rate of around 
0.1M E dd ([Revnivtsev et alJl200lHCornelisse et al.ll2003t 



iGallowav et al.l [2008D . Some degree of steady burning 
could also explain the scatter in the measured a values 
in our sample, but we cannot rule out the other mecha- 
nisms discussed here. 

6.3. Comparison with other bursters 

Finally, we investigate whether the accreted fuel con- 
finement proposed to explain the observed burst recur- 
rence times of 4U 1728-34, and other mechanisms dis- 
cussed here, may also be applicable to other bursters. 
In their study o f 9 bu rsters observed with BeppoSAX, 
Corncliss e et al.l (J2003) found that the burst rate in- 
creases slowly up to the peak rate, reaching it at approxi- 
mately 1 — 2xl0 37 ergs -1 (or ~ Q.05-0. 1M E dd for mixed 
fuel) , after which it sharply drops. Corncliss e et al.l 
( 2003) attributed such a behavior to the onset of stable 
hydrogen burning, which occurs at the transition from 
the lowest to the medium ac cretion rates (i.e., fro m the 
burst ignition regime 3 to 2; iFujimoto et al.lfl98~Th . 

The drop in the observed burst rate, at a slightly 
lower X-ray lumino sity, was also later confirmed by 
IGallowav et al.l (|2008l ). who studied a much larger sam- 
ple of bursters covering a much wider range of accretion 
rates. However. IGallowav et al.1 (J20081 ) argue that the de- 
creasing burst frequency is mor e likely due to onset of sta - 
ble burning of helium (see also lvan Paradiis et a.l.lll988T> . 
which must occu r at lower accret i on ra tes than theoreti- 
cally predicted. Gall oway et al.l ([2008D also found that 
the burst rate continues to decrease further until the 
bursting entirely stops (at ~ 0.3AfEdd)- In addition, 
the a parameter was found to remain constant or to 
slightly increase globally, rather than decrease as would 
be expected at the transition from the smallest accretion 
rate (ignition regime 3) to the medium accretion rate 
regime (2) and onset of stable burning of hydrogen (i.e., 
the transition of the H/He ignition to the He-ignition; 
IFuiimoto et~aTlll981[ ). 

Interpreting their results, iCornelisse et al.1 ()2003[ ) 
pointed out that there is a large discrepancy between 
the theoretically predicted (~ O.OlAfEdd) and mea- 
sured boundary between the two regimes, and offered 
an alternati v e exp lanation, which was proposed earlier 
by iBildstenl (120001) base d on t he llnogamov k. Sunvaevl 
(1999) model |Bildsten| (|2000D suggested that the ob- 
served burst rate drop is caused by the decrease of the 
local accretion rate. This would be possible if the nu- 
clear burning is restricted to a smaller area, which in- 
creases with the global accretion rate. If the area in- 
creases faster than the global accretion rate, the local ac- 
cretion rate may be decreasing. The correlation between 
the black-body emitting radius, measured in the burs t 
cooling tail of EXO 0748-676 ([Gottwald et all Il986h . 
and the corresponding accretion rate supports this hy- 
pothe sis. The observed a increase, following the burst 
drop | Galloway~ et al.l 12008). is also consistent with the 
Bildstcn (2000) interpretation (in this case the transition 
would be from the highest accretion rate regime (1) to 
regime 2, rather than from 3 to 2). Although the frac- 
tional covering of the accreted fuel may not be required to 
explain the bursting behavior of all sources, it certainly 
seems plausible at least for some bursters. 

Finally, we note that the same mechanism, which we 
believe is responsible for the large a scatter measured 
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for 4U 172 8-34 (i.e., the fuel sto rage and its ignition due 
to mixing; iFuiimoto et al.lH987J ). is also the most likely 
explanation for the large a scatter and weak frequent 
bursts with very short recurrence times of 10 minutes, ob- 
served in burste rs accreting mixed fuel (e.g., 4U 1636-53 , 
EXO 0748-676; ILewinet all [19871 IGottwald et aLllT986l 
1987). However, the scatter in a for the latter sources 
(10-500 for 4U 1636-53, and 5-50 for EXO 0748-676) is 
somewhat larger than that measured for 1728-34, likely 
because the fuel composition for sources accreting mixed 
H/He can vary due to steady H-burning prior to igni- 
tion. Interestingly, while our results suggest that some 
fraction of the burst fuel is typically left unburnt in 4U 
1728-34, short recurrence time bursts are not observed in 
this source, whereas the same mechanism is thought to 
give rise to such b ursts in systems accreting mixed H/He 
(jKeeket al.l[20T0h . 

7. CONCLUSIONS 

We compared the observed properties of 38 bursts de- 
tected in Chandra, RXTE and INTEGRAL observations 
of the helium-rich accretor 4U 1728—34 with new igni- 
tion models. We find that the observationaly-inferred ig- 
nition depths, assuming complete fuel spreading on the 
neutron star surface, are significantly smaller than the 
theoretically-derived minimum possible ignition depth of 
1 - 2 x 10 8 g cm 2 . 

One way to reconcile the observed and predicted burst 
recurrence times would be to assume that the observed 
X-ray luminosities underestimate the accretion rates 
(for example, due to a non-detection of the extremely 
soft or/and hard spectral components). However, for 



4U 1728—34 this scenario is not plausable because these 
additional spectral components would have to be 2—6 
times larger than the mesured bolometric luminosity. In 
addition, it would imply significantly larger a values than 
theoretically predicted. 

Alternatively, the ignition column could be increased, 
without increasing the global accretion rate, if we as- 
sume that the accreted material is confined to a rcstictcd 
area on the neutron star surface. To match the obser- 
vations with our ignition model, we find that the spot 
radii would have to be in the range 4.5—6 km, and that 
they seem to be cor related with t he glo bal accretion rate 
as predicted by the lLamb et al.l (|2009[ ) model. However, 
detailed ignition models that include all possible effects 
of the magnetic confinement are needed to confirm this 
result. An additional confirmation could come from the 
comparison of the X-ray emission from the accretion disk 
with the nuclear oscillation amplitudes detected simulta- 
neously from the source. 

To explain the weak, frequent bursts observed from 
4U 1728—34, we also consider shear-triggered mix- 
ing of the accreted helium to larger column depths 
(|Fuiimoto et al.l 119871: fPiro fc Bildstenl 120071 ). However, 
while this mechanism could push the fuel buffer from the 
outer envelope closer to the burning zone and explain 
the observed large a deviations, it might not be suffi- 
cient to explain the observed burst properties alone at 
the relatively low observed accretion rates. 

We would like to thank Ed Brown, Maurizio Falanga 
and Stratos Boutloukos for helpful discussions. We would 
also like to thank the anonymous referee for the useful 
and constructive comments. 
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Fig. 5. — Top: The inferred radius of the neutron star area cov- 
ered by the accreted material plotted versus the observed global 
accretion rate. The assumed neutron star radius is 10.4 km (Sec- 
tion 4) while the covered area was calculated to match the pre- 
dicted burst recurrence times. Bottom: The theoretical (full line) 
ignition column (calculated from the accretion rate; see also the 
middle panels of Fig. [4{ and the observed (la error bars) ignition 
columns (measured from the burst fluences; see the bottom pan- 
els of Fig. |4j plotted versus the scaled accretion rate. The scaling 
factors for each measurement are the same as in the left panel. 



